Selective and efficient quantum process tomography with single photons 
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We present the results of the first photonic implementation of a new method for quantum process 
tomography. The method (originally presented by A. Bendersky et al, Phys. Rev. Lett 100, 
190403 (2008)) enables the estimation of any element of the chi-matrix that characterizes a quantum 
process using resources that scale polynomially with the number of qubits. It is based on the idea 
of mapping the estimation of any chi-matrix element onto the average fidelity of a quantum channel 
and estimating the latter by sampling randomly over a special set of states called a 2-design. With a 
heralded single photon source we fully implement such algorithm and perform process tomography 
on a number of channels affecting the polarization qubit. The method is compared with other 
existing ones and its advantages are discussed. 
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The complete characterization of a linear quantum pro- 
cess mapping initial states into final states {£{pin) = 
pout) is a very hard task. In fact this well known fact 
can be understood as follows: A general completely posi- 
tive (CP) map can be fully characterized by the so-called 
X-matrix[l . Choosing the basis Ea consisting of op- 
erators, the x-matrix representation for the map S is such 
that S{p) = ^abXabEapEl {D = 2^ for a system of n 
qubits). The matrix Xab must be positive and hermi- 
tian for the channel to be completely positive (CP) and 
hermit ian itself. Moreover, for the map to be trace pre- 
serving its x~matrix should be such that the condition 
XabE^E}) = / is Satisfied. The matrix Xab is defined 
by D"^ — real parameters, a number that scales expo- 
nentially with the number of qubits of the system. This 
implies that full quantum process tomography (QPT) is 
unavoidably hard. Moreover, until recently, known meth- 
ods to estimate matrix elements Xab were also inefficient 
since they require resources (measured in terms of the 
number of repetitions of each experiment, on the number 
of operations required to perform them, etc) scaling ex- 
ponentially with the number of qubits [ln3| . Recently we 
introduced an efficient strategy enabling the estimation 
of any Xab element investing for such purpose resources 
scaling polynomially with the number of qubits of the 
system (and on the accuracy required for the estimation) 
[4]. In this paper we will present the first successful ex- 
perimental implementation of such strategy using single 
photons and linear optical elements. The experiment is 
not only an illustration of the use of a general scheme but 
also makes clear the advantages of the new tomographic 
method over its predecessors. 

It is convenient first to present the main ingredients 
of the only available method to perform selective and 
efficient quantum process tomography (SEQPT). For this 
purpose, it is best to first notice a general property of 
quantum channels. Assuming only that the operator base 
Ea is orthonormal in the Hilbert Schmidt inner product, 
it can be shown that any element Xab is related to the 



average fidelity of a modified channel Sab- Such channel 
depends on the original S and on the operators Ea and 
El). It acts on any state p as £ab{p) = £{EapEi)). Thus, 
we can show that [4j 

Xab = ^iiD+ l)FiSab) - Sab) ■ (1) 

Here, the average fidelity F{£ab) is computed integrating 
over all pure states using the (unitarily invariant) 
Haar measure, i.e. 

F{£ab) = J d\<P){<P\£ab{\4>){m4>)- (2) 

The method of SEQPT is based on the use of the above 
identity and on the fact that the integral over the en- 
tire Hilbert space can be exactly computed using finite 
resources. It has been shown that the average fidelity 
can be exactly evaluated by computing the average over 
a finite set of states that form a so-called 2-design [5] . In 
fact, if the set S = {IV^j), j = 1, •.•,^} is a 2-design, we 
have 

Fi^-b) = ^ Y.('^j\£ab{\^j){mi'j)- (3) 

j 

Although the exact computation of Xab using a 2-design 
would involve performing a sum over an exponentially 
large set, the estimation of such coefficients with finite 
(D-independent) precision is possible by randomly sam- 
pling over states As described in [4 and discussed 
with more detail below, the precision in the estimation 
scales better than 1 / ^/M with the number of repetitions 
M of the experiment. 

There are several known examples of 2-designs such as 
the one formed by the D{D-\-l) states belonging to 1^ + 1 
mutually unbiased bases [6 (two bases are mutually un- 
biased iff the absolute value of the overlap between states 
in such bases is 1/^/D). The equation (3) suggests an im- 
mediate way to characterize the quantum channel: For 



every coefficient Xab one needs to compute the survival 
probability of states when evolved over the channel 
Sab and then average over the 2-design. For the strategy 
to be viable one needs not only to make sure that the 
channel can be implemented efficiently but also that the 
sum over the elements of the 2-design can be estimated 
with polynomial resources. The fact that these two tasks 
can be performed is the core of SEQPT. In particular, 
the diagonal elements Xaa are average fidelities of the 
CP-map £aa^ obtained by composing the original map 
S with the unitary operation Ea'. £aa{p) = £{EapEa). 
Evaluating off diagonal elements Xab requires a different 
strategy that at first sight seems to be rather different 
than the one used for diagonal coefficients. However, the 
main point in this strategy, again, is to realize that any 
off diagonal coefficient Xab is related through (1) with 
the average fidelity of the map Sab- This map is not CP 
but can be obtained as the difference between two CP 
maps. For this reason, as shown in [4], the off diago- 
nal coefficients can be obtained as the mean value of an 
ancillary qubit conditioned to the survival of the state 
l^pj) and averaged over the 2-design. The complete al- 
gorithm is described in detail in Figure 1. In fact, real 
and imaginary parts of Xab are obtained by conditionally 
measuring the mean values of ax and ay of an ancillary 
qubit that interacts with the system with controlled-£^a^5 
operations. 
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FIG. 1: The quantum algorithm for measuring Re(xab) for 
a given channel 8. The method requires an extra ancillary 
(clean) qubit. The imaginary part of Xab is estimated in the 
same way by measuring the polarization of the ancilla along 
the y axis. 

Although SEQPT is the only available scheme to ef- 
ficiently estimate any coefficient of the chi-matrix it has 
never been implemented in practice. Other methods have 
been demostrated in different experimental setups (see, 
for example, [7HT2]). We present here the results of the 
first implementation of SEQPT in a photonic quantum 
information processor. For this we use a heralded single 
photon source encoding two qubits in a single photon. We 
do this by using both the polarization degree of freedom 
and the momentum (path) degree of freedom of a her- 
alded photon generated by parametric-down-conversion 
[13]. The method enables us to characterize unknown 
quantum channels affecting the polarization qubit using 
the path qubit as an ancilla. 

The full implementation is schematically shown in Fig- 
ure [2] For the sake of clarity we divide the description of 
the method in four steps: i) generation of initial states of 
the 2-design ii) application of controlled operations 
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FIG. 2: A continuous laser diode at 405nm and 40mW gen- 
erates frequency degenerate twin photons on BBO crystal, 
cut for Type-II parametric-downconversion. One photon is 
used as a herald while on the other a polarization and path 
qubit are encoded. We first prepare a state from the 2- 
design on the polarization qubit with a polarizer, half and 
quarter wave-plates, then we perform a Hadamard gate on 
the path with a non-polarizing 50:50 beam-splitter, next we 
perform the controlled operations Ea and Eb with a quarter- 
half-quarter wave plate sequence on each arm, both paths 
are then sent through the unknown channel S affecting the 
polarization qubit equally on both paths. Finally measure- 
ment of ax or ay in the path qubit is done by interfering both 
paths at a second non-polarizing beam-splitter with a phase 
previously set by tilting a glass; projection on the prepared 
polarized state is done by a quarter and half wave-plate fol- 
lowed by a polarizing beam-splitter and detection is done by 
fiber coupled avalanche photodiodes. The interferometer has 
an average visibility of 92%. Not shown in the figure is an in- 
tense beam counter-propagating in the interferometer which 
is monitored to actively control the position of one of its mir- 
rors with a piezoelectric disc attaining a stability of over A/30 
over all measurements. 



between ancillary qubit and target qubit , iii) the evolu- 
tion through the unknown process and iv) measurement 
in the ax or ay base for the ancilla (path) qubit con- 
ditioned on the survival of the (polarization) states. 
With these steps we implement the algorithm of Figure [l] 
is with the photonic setup described in Figure [2] Stage (i) 
is rather simple. For the case of a single qubit it turns out 
that the 2-design is formed by the D{D -\- 1) = 6 eigen- 
states of the three Pauli operators (which define three 
mutually unbiased bases). For the polarization of pho- 
tons these are: two states with vertical-horizontal polar- 
ization, two with diagonal (45-135 degrees) polarization 
and two with circular (right-left) polarization. The gen- 
eration of such states is done with the standard polarizer, 
half and quarter wave plate configuration. Stage (ii) of 
the algorithm is conceptually simple: the photon is split 
by the first non-polarizing beam-splitter, which acts as 
the Hadamard gate in the ancillary qubit, then any con- 
trolled operation on the polarization qubit can done by 
rotating independently the polarization in each path us- 
ing a quarter-half-quarter wave plate configuration [13]. 
In order to apply the unknown process (iii) both paths 
are then made parallel and sent though a zone where 
the unknown channel £ is performed on the polarization 
qubit. Finally (iv) to measure the expectation values of 
ax and ay conditioned to the survival of we make 
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both paths interfere with a relative phase ^ at a sec- 
ond non-polarizing beam-splitter by tilting a glass in one 
of the paths and project into the polarization states with 
the inverse quarter, half wave plate and polarizer configu- 
ration as in (i). With this scheme measure the quantities 
Pab{o'x{y)^^'->^tl^j)^ which are the probabilities of finding 
the ancilla in the ± state of the x{y) base conditioned on 
the survival of the input states {i/jj) for each Sab- With 
these probabilities we obtain all the necessary data to 
determine the matrix element Xab- A note must be made 
about diagonal elements: as mentioned above, in this 
case the scheme simplifies significantly because there is 
no need for an ancilla qubit so one then can look at only 
one arm of the interferometer by blocking the other one. 

The elements of Xab were detected for two different 
noisy processes affecting the polarization degree of free- 
dom. The results are shown in Figure [3] We must stress 
that although we measured all the matrix elements char- 
acterizing both channels our method determines any Xab 
independently and efficiently (i.e., as opposed to previous 
methods [l] it is not necessary to fully characterize the 
channel in order to determine a single Xab coefficient). 
The measured processes correspond to the identity chan- 
nel (i.e., free propagation through air) and to a quarter 
wave plate at 0°. Those processes were also fully char- 
acterized by means of the standard method of QPT (as 
explained in chapter 10 of [1 ). The results obtained by 
both methods are compared in Figure [3] and turn out to 
be in very good agreement. As a figure of merit to com- 
pare both schemes we numerically calculate fidelity be- 
tween the channels determined by each method obtaining 
F = 95, 1%±1.5% for the identity and F = 96, 3%±1.6% 
for the half wave-plate at 0° [16]. The maximum fidelity 
attainable is mainly limited by interferometer visibility, 
which results in a not perfectly clean ancilla. 

SEQPT is a method that is suited to perform par- 
tial process tomography selecting the relevant param- 
eters one wants to estimate and investing polynomial 
resources for such estimation. It is interesting to com- 
pare the resources required to implement this new tomo- 
graphic method and previously existing ones. To deter- 
mine any single matrix element Xab using SEQPT, we 
measured D{D -\- 1) = 6 survival probabilities. On the 
contrary, obtaining a single matrix element Xab using the 
standard QPT [1] requires estimating D^xD^ = 16 tran- 
sition probabilities. This implies that even at the level 
of a single qubit the SEQPT is more "efficient". This 
comparison might seem somewhat unfair since after such 
number of experiment the standard method provides all 
the information required to estimate the full Xab matrix 
and not just a single element of such matrix. However, 
when working with larger systems, it is clear that one 
would never estimate the exact matrix elements Xab by 
performing the average over the entire 2-design. In fact, 
this would always be exponentially hard (as the 2-design 
is a set containing an exponentially large number of el- 
ements). Instead, the main virtue of SEQPT is that it 




FIG. 3: Experimental results for SEQPT (a.l) and (a.2) and 
Standard QPT (b.l) and (b.2). First two rows display mea- 
sured real (top) and imaginary (bottom) parts of the matrix 
Xab of the quantum process corresponding to free propaga- 
tion (identity channel). Similarly the last two rows display 
correspond to real (top) and imaginary (bottom) parts of the 
matrix Xab measured for a process corresponding to a QWP 
at 0°. 

enables us to estimate any element with fixed accuracy by 
performing a number of measurements that only depends 
on such accuracy and is independent on the number of 
qubits. As shown in [4] the error in the estimation of 
the average in equation (3) after M experiments scales 
as 1/a/M. However, it is interesting to get a tighter 
bound for such error. In fact, if the average fidelity is 
computed by sampling over states {i/jj) which are never 
prepared twice, then one can show that the error in the 

average turns out to scale as Ax cx ^-^(1 — ) where 
K = D{D-\-l) is the number of elements in the 2-design. 
In fact, this equation tells us that the error vanishes as 
the number of experiments M = K. 

Having the experimental data at hand we can test the 
behavior of the estimation error. Using the raw data 
we split the 2-design in samples of different (variable) 
size and computed the average over each of such sam- 
ples. The results are shown in Figure |4] where we display 
the behavior of the estimation error for the matrix ele- 
ment Xab obtained by sampling the 2-design in groups of 
increasing number of states (the size of the sample, M 
grows up to K = D{D -\- 1) = 6, which is the cardinal of 
the 2-design). These results are quite interesting and not 
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FIG. 4: Convergence for the chi matrix elements theoretical 
(smooth line) and experimental (segmented lines). Shown are 
the errors in the estimation of every Xab for both measured 
processes. Each curve shows the error made in the estimation 
for each Xab if only M measurements were made. All possible 
choices of the first M states are plotted and they lie below the 
theoretical bound. The insets show the errors for two specific 
Xab elements. In (a) the error can stay close to the bound 
for some choices of the first M measurements while in (b) is 
always almost null. 

expected: Although for a single qubit the sample space is 
rather small the above mentioned bound turns out to be 
satisfied rigorously for every instance of the experiment 
and not only at the statistical level. We see from Figure 
[4] that all possible errors on all possible choices of sample 
partitions lie below the analytical bound. It is clear that 
this is not expected for all possible random distributions, 
but the ones realized in the experiment strictly satisfy 
the bound. In fact, it is not hard to imagine possible 
values for the results of the experiment that would vio- 
late the bound for certain samples but hold to it at the 
statistical level. However, such cases are not realized in 
the experiment which suggests that it would be possible 
to find a tighter bound. 

We conclude by noticing that looking at the quantum 
state of the photon at the output of the final beam splitter 
in Figure 2 is quite revealing about the nature of SEQPT. 
In fact, it can be shown that if we prepare the initial state 
of the polarization qubit to be then the unnormal- 
ized output state turns out to be described by the density 
matrix £{{Ea ± Eij)\\iljj){iljj\{Ea ± £^5)) (the ± signs are 
controlled by a = 0, tt phase with the glass in Figure 2). 
This shows that the ancillary (path) qubit plays a rather 
simple role in SEQPT: it is simply a tool to prepare a spe- 
cial initial state which is later sent through the channel 
£. This observation serves as a motivation for a simple, 
but illuminating, reformulation of the SEQPT algorithm. 
In fact, SEQPT is connected with a class of quantum al- 
gorithms known as "deterministic quantum computation 
with one clean qubit" (DQCl) [T4l|T5]. However, as op- 
posed to DQCl, Figure 1 shows that the measurement 



on the ancilla can be moved to the initial stage of the 
algorithm (i.e., it can be performed before the channel £ 
acts on the system). From this point of view, SEQPT 
can be described as follows: The unnormalized state of 
the system conditioned on the detection of the ±1 eigen- 
value in the measurement of ax on the ancilla turns out 
to be {Ea ± Ej,)\i;j){i;j\{Ea ± E^). Thus, after the sys- 
tem goes through the channel f , the probability of de- 
tecting the state {t/jj) is simply the fidelity of the map 
£±{p) — £{{Ea i Eij)p{Ea ± ^5)). The expectation value 
of ax is computed as the difference between the fidelities 
of the maps £± , which is nothing but the real part of Xab 
(analogously for the imaginary parts with (j) = ±7r/2). 
This new interpretation of SEQPT, suggested by the ex- 
periment indicates that the method could indeed be im- 
plemented without ancillary qubits at all. In any case, 
the present experiment shows the viability of SEQPT as 
a way to characterize a channel affecting polarization of 
single photons and displays its advantages over other al- 
ternative approaches. 

The experiments were performed in the quantum op- 
tics lab at CITEDEF. JPP and MAL are members of 
CONICET. The authors acknowledge discussions with 
A. Bendersky, A. Hnilo and R. Piegaia. 
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